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by 
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ABSTRACT* 


The  problem  studied  Is  that  of  a  slab,  heated  in  an  arbitrary  manner 
on  one  face  and  insulated  on  the  other^  which  melts  (or  solidifies )>  the  loate- 
rial  being  allowed  to  restain  stationary  after  change  of  phase.  Variable  ma¬ 
terial  properties  are  taken  into  account.  After  preliminary  general  consider- 
atlonSj  it  is  shown  that  the  solution  to  the  stated  problem  is  unique.  It  is 
then  proved  that  higher  rates  of  melting  and  higher  tenperatures  will  result 
from  certain  combinations  of  the  magnitude  of  the  applied  heat  input  and  of 
a  fictitious  heat  source  traveling  with  the  solid-liquid  interface.  From  this 
result  a  method  is  developed  for  the  construction  of  upper  and  lower  bounds  to 
the  solution  of  the  problem;  an  example  is  also  presented.  It  is  also  shown 
that,  under  the  same  arbitrary  heat  input,  the  rate  of  melting  in  the  present 
problem  is  always  lower  than  that  in  the  companion  problem  in  which  the  mate¬ 
rial  Is  instantaneously  reoioved  after  change  of  phase. 


*  This  work  was  performed  as  part  of  a  project  sponsored  by  the  Office  of 
Naval  Research. 
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1.  Introduction 

The  solution  of  heat  conduction  problems  of  fltlng  or  solidification 
often  presents  considerable  mathematical  difficulties,  and  Is  therefore  often 
approached  cither  by  numerical  or  approximate  techniques.  A  different  approach 
was  described  In  [ll,  in  which  a  method  of  constructing  upper  and  lower  bounds 
for  the  temperature  and  for  the  unknown  location  of  the  melting  front  In  this 
typo  of  problem  was  devised.  The  problem  considered  there  was  that  of  a  slab, 
Initially  solld^.  Insulated  on  one  face  and  subjected  to  an  arbitrary  heat  Input 
on  the  other;  after  melting  starts,  the  known  heat  Input  Is  applied  directly 
to  the  moving  boundary  of  the  solid.  In  the  present  work  the  same  type  of  ap¬ 
proach  Is  extended  to  Include  the  problem  of  a  melting  slab.  In  which  the  melted 
material  Is  not  removed,  but  remains  stationary. 

Many  of  the  theorems  proved  In  [1]  will  be  needed  In  the  course  of  the 
present  proofs,  and  are  therefore  restated  here  (with  some  minor  extensions)  for 
the  sake  of  convenience,  at  the  beginning  of  Section  2.  The  principal  results 
required  for  the  establishment  of  bounds  In  the  present  problem  follow  In  the 
same  section,  while  a  statement  of  this  particular  melting  problem  and  a  proof 
of  uniqueness  of  Its  solution  are  given  In  Section  5*  Two  types  of  upper  and 
lower  bounds  are  established  In  Section  If:  the  first  compares  melting  rates  and 
tenqieratures  in  the  present  problem  under  certain  combinations  of  the  magnitude 
of  the  applied  heat  Input  and  of  a  fictitious  variable  heat  source  traveling 
with  the  liquid-solid  Interface,  while  the  second  compares  the  solution  of  the 
present  problem  with  that  of  [1].  The  construction  of  the  first  type  of  these 
bounds  and  their  use  In  obtaining  estimates  of  solution  Is  discussed  together 
with  an  Illustrative  example.  In  Section 

^  The  work  of  [l],  as  well  as  the  present  one,  apply  equally  to  the  problems 
of  melting  and  of  solidification;  for  convenience,  however,  only  Che  former 
type  of  problem  will  henceforth  be  referred  Co. 
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2.  Baalc  Theorems 

The  first  two  of  the  following  theoreos  are  listed  here  for  ease  of  ref¬ 
erence;  their  proofs  laay  be  found  in  [1]. 

(2) 

Theorem  I  Let  u(p)  be  a  solution,  in  class  C'  of  the  equation 

a(P)^  +  b(p)^  -  1^  -  0  ;  a  i  0  (l) 

dx 

at  all  points  F(x>t)  of  a  domain  D  in  the  xt -plana  where  a  and  b  are  real  fi¬ 
nite  continuous  functions,  and  let 

^  .  f(P),  P  on  {2a) 

X 

I 

u  ■  k  at  some  point  on  B-B_^  (2b) 

where  n^  denotes  the  direction  of  the  component  parallel  to  the  x-axls  of  the 
interior  normal  to  the  boundary  B,  k  is  a  constant, and  where  B_^  is  that  part 
of  B  which  includes  all  points  P  such  that  (a)  the  interior  normal  to  the 
boundary  exists  and  Is  directed  In  the  negative  t  direction,  and  (b)  each 
point  P  is  an  interior  point  of  B_^.  Then: 

(a)  if  f(P)  ■  0,  u  -  k  throughout  D; 

(b)  if  f(P)  is  prescribed  throughout  B-B_^,  then  u  Is  uniquely  deter¬ 
mined  throughout  D. 

2 

Theorem  II  Consider  a  simply  connected  domain  D  In  the  zt-plane, 
whose  boundary  is  (Fig.  l)  formed  by  segments  of  the  straight  lines  x  ■  x^,  and 

p 

This  Is  a  slight  extension  of  Theorem  II  of  [1]. 


-k- 


t  ■  Cg  and  by  a  line  defined  by  a  continuous  single'^alued  function  x  »  F(t) 
satisfying  Llpschltz  conditions  for  t^  <  t  <  t^  and  intersecting  the  line 
t  s  tj^  at  a  point  :  (x^^t^).  Let  u(P)  be  the  solution  of  Eq.  (2)  In  D  with 


■  0  on  X  a  X. 
ox  0 


U  =»  “q(*)  t  ■  tj^ 


(3a) 


and 


either  «  -f(t)  or  u  ■  f(t)  on  x  ■  F(t) 


(5b) 


Then,  in  either  of  these  cases,  if  f(t)  i  0  and  Uq(x)  ^  0,  u  2  0  throughout  D, 
and  If  f(t)  ^  0,  and  Uq(x)  ^  0,  then  u  s  0  throughout  D. 

Theorem  III  Let  Uj^(P)  and  U2(P)  be  solutions  of  Eq.  (l),  respectively 
in  the  domains  Dj^  and  Og  of  the  xt*plane  defined  in  Fig.  2,  and  let 


du.(x,,t) 

sr  ■  '  'o  ‘ h 


(4a) 


^l(xi+L,t) 

^  1  .0,  tQ<t<t2 


(4b) 


dT  -  ^  h  <  <  *^2 

'»l(x,to)  -  0  ,  Xj^  <  X  <  L 


(4c) 

(4d) 


du. 


M 


Ui  -  Ug 


dx 


on  X 


F(t)  ,  t,  <  t  <  t. 


(4f) 
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where  M(P)  >  0  end  N(P)  >  0  and  where  F(t)  Is  a  given  Llpschita  continuous 
single-valued  function  of  t*  Then 

(a)  if  f(t)  ■  0  in  t^  <  t  <  tg,  Uj^  ■  0  and  Ug  a  0  throughout  their 
domains  of  definition; 

(b)  If  f(t)  Is  arbitrarily  prescribed  In  t^  <  t  <  tg,  then  u^  and  Ug 
are  uniquely  detertalned  throughout  their  respective  domains; 

(c)  if  f(t)  *  0  in  tQ  <  t  <  tg,  then  Uj^  *  0  and  Ug  S  0  in  their  re¬ 
spective  domains,  and 

(d)  if  f{t)  i  0  in  t^  <  t  <  tg,  then  Uj^  i  0  and  Ug  s  0  in  their  re¬ 
spective  domains  of  definition. 

Statement  (b)  is  a  corollary  of  statesient  ^a),  and  it  %d.ll  be  apparent 
that  the  proof  of  (c)  requires  only  obvious  otodificationa  to  apply  to  part  (d) 
as  well.  Hence  only  parts  (a)  and  (c)  need  be  considered  below. 

Proof  of  (a)  Clearly  (Theorem  I)  u^^  ■  0  for  t^  <  t  <  t^^.  For 
t,  <  t  <  t„,  ^1  and  ^2  have  opposite  signs  or  vanish  on  a  >  F(t).  If 

they  are  both  zero.  Theorem  I  again  applies  and  Uj^  ■  0,  Ug  ■  0.  If  they  were 
not  identically  zero,  then  on  x  -  F(t)  one  cculd  find  values  t'(2  t^)  and 
ft  >  0  such  that 


tj^  <  t  <  t' 


and 


either 


or 


0  , 

0  , 


_2 

dn 


X 


X 


<  0 


>  0 


t'  <  t  <  t*  +  6 


(5a) 


(5b) 
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Tben,  by  Theorem  II,  for  the  first  of  these  possibilities  ^  0  and  ^  0  in 

t'  <  t  <  t'  +  6  ;  therefore,  on  x  -  F(t),  Uj^  *  ■  0  in  this  interval  of  time, 

and  similarly  for  the  second  of  the  possibilities  in  Eqs .  (^).  But  then 

Theorem  II  [applied  first  to  the  domain  t'<t<t'  +  6,  Xj^<x<  F(t)  and  then 

to  the  domain  t'  <  t  <  t'  +  B,  F(t)  <  x  <  Xj^  +  1.]  shows  that  in  either  case 

a  3  0  at  any  point  within  this  time  interval,  so  that  neither  of  (^b)  can 

exist.  Hence  part  (a)  of  the  theorem  is  proved. 

Proof  of  (c)  From  Theorem  II,  Uj^  ^  0  in  t^  <  t  ^  tj^;  for  tj^  <  t  <  tg 

two  possibilities  can  arise,  depending  on  whether  initially  ^1  is  negative  or 

dx 

positive  on  X  a  F(t)  in  this  interval.  In  the  former  case  Eqs.  (4f)  give 

i:  0  on  X  ■  F(t)  and  hence  Ug  ^  0  in  domain  2,  and  thus  <  0  on  x  >  F(t). 
Hence,  by  Theorem  II,  u^  is  non-positive  throughout,  as  was  to  be  proved.  Sup¬ 
pose  now  that  a  later  time  t*  exists,  at  which  ^1  on  x  ■  F(t)  changes  sign, 

dx 

or,  in  other  words : 


Ui  = 

Ug  ^  0 

_ 1 

dx 

dx 

du. 

X 

dx 

>  0  j 

tj^  <  t  i  t' 


t  -  t' 


throughout 


du. 


(6) 

on  X  ■  F(t) 


0  •  >  0  t'<t<t*+6.  6>0 

'ox  ' 


Then  (Theorem  II)  Uj^^0  1nt'<t<t'  +  B  for  all  x  and  so  Ug  SO  in  this 
time  interval  on  x  >  F(t).  Now  if  Ug  were  non-positive  on  x  ■  x^,  Theorem  I 
would  insure  that  it  be  non-positive  throughout,  as  was  to  be  proved;  if,  on 
the  other  hand,  Ug  were  positive  on  x  «  Xj^  at  some  time  t"  (t'  <  t"  <  t*  +  6), 
then  some  time  within  t'  <  t  <  t"  the  condition  ^2  >  0  as  well  as  (to  avoid 
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an  interloc  ii«ucli»im  m  required  by  Plcoae's  theorem  [l»2l)  the  coadltlone 
■  0  and  — ^  i  0,  all  on  x  ■  Xp  would  have  to  prevail*  The  conditions 


are  however  in  contradiction  with  Eq.  (l),  and  therefore  nust  be  non-poei- 

tive  again.  Now  it  may  happen  that  at  a  later  time  ^1  on  x  «  F(t)  once 

dx 

more  changes  sign:  then  the  first  part  of  this  proof  applies  once  again,  and 
similarly  for  any  further  changes  in  sign.  There  now  remains  the  second  pos¬ 
sibility  alluded  to  at  the  beginning  of  this  proof,  namely  that  ^1  is  ini- 

dx 

tially  positive  on  x  »  F(t);  in  this  case  the  second  part  of  the  above  proof 
applies  first,  and  any  later  changes  in  sign  of  this  derivative  can  be  taken 
care  of  as  above.  Then  proof  is  thus  coopleted. 


3*  Statement  of  the  Melting  Problem^ 

Consider  a  slab,  initially  (i.e.  at  t  >  o)  solid  at  zero  temperature 

and  occupying  the  region  0  <  x  <  L,  and  iimulated  at  x  ■  1>.  An  arbitrarily 

prescribed  heat  input  Q(t)  is  applied  at  x  ■  0,  so  that  the  temperature  in  the 

slab  rises  and  at  x  ■  0  reaches  the  melting  temperature  T  at  the  time  tut. 

n  m 

Melting  continues  to  take  place  thereafter,  and  a  portion  of  thickness  8(t) 
is  taken  to  have  melted  at  any  time  t  2  t^,  while  the  prescribed  heat  input 
Q(t)  still  continues  to  be  applied  at  x  ■  0.  Thus  at  any  time  t  >  t^,  the 
portion  of  the  slab  within  s(t)  <  x  <  L  is  still  solid,  while  that  within 
0  <  X  <  s(t)  is  liquid;  the  subscripts  S  and  L  will  be  used  in  what  follows  to 
distinguish  quantities  pertaining  to  the  two  phases  of  the  material.  The  math¬ 
ematical  fomulation  of  the  problem  is  as  follows  for  the  temperature 

^  A  statement  and  discussion  of  this  and  of  the  analogous  solidification  prob¬ 
lem  were  given  for  example  in  [3],  An  approximate  solution  of  this  problem 
was  presented  in  18]. 
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T-(xft)  <  T  In  the  solid: 
S  B 


^(k  ^)- 

ax  ax  ^ 

«  s 

Ps®s  at  ' 

e(t)  <x<L,  0<t<tj^ 

(7) 

Tg(x,0)  -  0 

0  <  X  <  L 

(8) 

3Tc 

-rf(L,t)  - 

0  i 

0  <  t  <  tj^ 

(9) 

■*‘s  ax 

-  Q(t)  ; 

0  <  t  <  t 

m 

(10) 

end  for  the  temperature  T, (x,t)  k  T  in  the  liquid: 

JL  ul 


ar 


ar 


ar 


(11) 

(12) 


with  the  following  interface  conditions: 


Tg[8(t),tl  -  Tj8(t),tl  - 

V  ^(8(t),t]  5s[.(t),t] 

Tan  dx  Sm  dx 


t  <  t  <  t-  /,_x 
m  L  (13) 


«•(')  it‘'’ 


and  with 

8(t)  .  0  ;  0  s  t  s  t^  ilk) 

The  tines  t  and  t,  are  respectively  defined  by  the  equations 

Ql  la 
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Kid  s(tj^)  -  L  (15) 

The  thermal  dlffualvlty  h  ■  k/(pc),  the  conductivity  k,  the  specific  heat  c, 
and  the  density  p  are  asauaed  to  be  functions  of  the  temperature  and  therefore 
vary  with  both  z  and  t.  The  subscript  m,  affixed  to  any  quantity,  indicates 
that  the  value  at  the  melting  temperature  T^  must  be  used.  The  latent  heat 
of  melting  is  denoted  by  I, 

The  heat  input  Q*(t)  appearing  in  Eqs.  (I3),  has  the  physical  meaning 

of  a  variable  heat  source  traveling  with  the  interface,  and  is  identically 

zero  in  the  problem  whose  solution  it  is  desired  to  find;  nevertheless  it  will 

be  included  in  all  the  derivations  because  it  is  a  convenient  quantity  to  deal 

with  in  the  calculation  of  bounds  for  the  solution  of  this  problem.  Q*  has 

not  been  defined  for  t  <  t  ;  it  will  be  convenient  to  take  it  as  zero  in  this 

m 

range. 

k 

It  can  be  readily  shown  by  means  of  an  overall  beat  balance  that 
t  s(t)  L 

/  [Q(t)  +  Q*(t)]  dt  -  /  ILdx  +  /  H  dz  +  (p_  i  +  H  )  s  (l6) 

0  0  ^  s(t)  ®  Sm  m 

where  the  heat  contents  and  are  defined  as 

-  /  Pg(T')cg(T')dr’  (17a) 


The  derivation  is  very  similar  to  that  of  [4],  [3I  [61;  the  extension  to 

the  present  case  of  variable  properties  presents  no  difficulty  since  in  terms 
of  the  heat  contents  and  the  right-hand  sides  of  Eqs.  (7)  and  (ll)  re¬ 
duce  respectively  to  (dH^/dt)  and  (dH^^/dc). 
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and 

T 

-  /  P,(T')cjT')dt'  (17b) 

T 

B 

Note  that  the  integrands  of  Eqs.  (I7)  are  positive  and  therefore  both  and 
are  monotonically  increasing  functions  of  T(x,t).  The  syobol  stands  for 

"sftj- 

Theorem  IV i  Uniqueness  of  Solution  It  will  now  be  proved  that  there 
exists  at  most  one  solution'^  to  Eqs*  (7)  to  (14),  corresponding  to  prescribed 
functions  Q(t)  and  Q*(t).  To  prove  this,  assume  that  two  distinct  solutions 
exist,  and  denote  them  by  the  subscripts  1  and  2;  then  Theorem  III  insures 
uniqueness  if  s^  ■  a^.  Only  the  possibility  s^^  ft  Sg  need  therefore  be  con¬ 
sidered,  or,  without  loss  of  generality,  we  may  set 


*l“*2*  O^tst' 

Sg>ej;,  t'<tit'  +  6,  6>0 


(18) 


Since  the  solution  is  known  to  be  unique  before  melting  starts  [4,71  b  is  the 

B 

same  for  each  solution;  hence  t'  2  t^.  After  the  start  of  melting,  write  Eq. 
(16)  for  each  solution  (at  some  time  t"  in  the  interval  where  Sg  >  s^)  and  sub¬ 
tract  the  results  to  get 


•l(t") 

; 

0 


L 


S2"*^S1^^ 


+  ^<»Sb^  +  Hj[s2(t")-sj^(f)l  (18a) 


a  . . . 

^  That  is,  a  twice  continuously  differentiable  function  T(x,t)  and  a  Lipschits 
continuous  function  s(t)*  Except  in  the  special  case  of  Eqs*  (24),  this  char¬ 
acter  of  the  solution  is  assuaed  throughout  the  reisainder  of  this  paper. 
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On  X  -  and  a  while  on  x  -  SgCt),  Tj^  ^  and  -  T^; 

hence  Theorem  111^  applied  In  turn  to  the  regions  Sg  <  x  <  L  and  0  <  x  < 

Inaures  that  3:  in  these  regions,  so  that  the  first  and  third  integrals 

in  (I8a)  are  non-negative.  Furthermore 

/hj^dx  i  0  and  /h^^dx  s  [s2(t")-Sj(t'')]  (I8b) 

»1  *1 

Hence  the  right-hand  aide  of  (l8a)  cannot  be  zero,  this  equation  cannot  be 
satisfied,  and  uniqueness  is  assured. 

U.  tinner  and  Lower  Bounds 

Theorem  V  Consider  two  solutions  of  Eqs.  (7)  to  (14),  denoted  by  the 
subscripts  1  and  2,  corresponding  respectively  to  heat  inputs  Qj^(t),  Q^(t)  and 

»  Q*  ;  QgCt)  +  QgCt)  i  Qi(t)  +  Q^(t)  (19) 

Then 

S2(t)  i  Sj^(t)  ;  t  2  0  (19a) 

T2(x,t)  i  T^(x,t)  ;  OixSL  (19b) 

It  is  understood,  in  Eq.  (19b)  and  elsewhere,  that  for  the  temperature  T  one 
must  read  the  pertinent  one  of  the  functions  T^  and  T^. 

The  case  in  which  in  both  Eqs.  (19)  the  equal  signs  hold  thrmighout  is 
covered  by  Theorem  IV  aid  need  not  be  considered  further.  To  prove  the  present 
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Cheorem  In  Che  general  case.  It  may  first  be  assumed^  without  loss  of  general- 
tty,  Chat  a  time  t*  exists  such  that 


QjCt)  +  Q*(t)  -  Qg(t)  +  Q*(t)  for  0<t«t' 

Ql(t)  +  Qj(t)  <Q2(t)  +  Q*(t)  for  t'  <  t  S  t'  +  6  ,  6  >  0 
It  will  now  be  shown  Chat  a  tine  t"(t'  <  t"  <  t*  -t-  6)  exists  such  that 
Sg(t)  >  Sj(t)  for  t'  <  t  S  t" 


(20) 


(20a) 


Assume  in  fact  that  this  la  not  so;  then  the  situation  Is  that  shewn  In  Fig. 

3a  (note  that  It  Is  Insaterlal  In  Che  proof  whether  t'  >  t^  or  t'  <  C^^^  though 
only  Che  former  case  Is  shown  In  Che  figure).  Eq.  (16)^  written  for  each  so¬ 
lution  at  t  ■  c",  gives  after  subtraction^ 

t"  ®2(^")  *l(^") 

0  <  /  [(Qg  Qp-(Qi  +  Q*)ldt  -  ;  (H^-H^i)dx  W^^^^^(H32-HLl)d* 

*  .A«S2’«Si)'^  -  (Psn^  *  H^)[si(f)-8  (f)l  (21) 

•  l(t  ) 

Howeveri  on  x  -  s^,  ■  ^n,  '^2  *  ^m*  while  on  x  ■  Sg,  Tg  m  and  Tj^  * 

hence  the  first  and  third  Integrals  on  the  right-hand  side  of  (21)  are  negative. 
Furthermore, 


■l  “l 

/  Hj^jdx  »  0  and  /  Hggdx  *  Hjs^  -  Sg) 

*2  *2 


(21a) 


-1^- 


Hence  the  right-hend  aide  of  (21)  cannot  be  positive  as  required^  and  (20a)  Is 
proved.  Note  that  the  validity  of  Inequality  (l9b)  for  the  range  of  times  con¬ 
sidered  has  also  been  proved. 

It  has  been  thus  proved  that  Initially  s^  cannot  exceed  Sg^  and  It  oust 
be  sho%in  next  that  no  later  time  can  exist  at  which  s^^  exceeds  Sg*  Indeed,  If 
there  were  such  a  time,  there  would  have  to  be  at  least  one  Instant  at  which 
Sg  ■  s^;  let  the  first  of  these  Instants  be  t'"  (>  t'),  as  shown  In  Fig.  3b. 

Note  that  then 

T.  -  T  -  T  ;  s,  -  s.  ;  i  ffi  at  P»  (22) 

^  dt  dt 

The  second  of  Eqs.  (13),  written  for  each  solution  at  t**,  gives,  after  subtrac¬ 
tion 

where  dots  Indicate  differentiation  with  respect  to  time.  Now  ^  0 

on  X  >  s^  and  hence  also  to  the  left  of  P^;  In  view  of  (22),  then, 

4^^L2'^L1>  ^  °  **=  ^2  ^22b) 

Similarly,  ‘^32'^Sl  ^  ^  ^  *  *2  hence  also  to  the  right  of  P^,  so  that 

^Vsl)  *  0  *“2 

This  means  that  the  left-hand  side  of  (22a)  Is  non-positive;  however  the  right- 
hand  side  of  this  equation  Is  noo-negatlve  In  view  of  (I9)  and  (22)  and  there- 
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fore  this  implies  a  contradiction  If  at  least  one  of  the  Inequalities  s^  < 
or  0^  >  Q*  hold.  In  the  special  case 


at  t  -  t'" 


(25a) 


the  above  proof  falls  since  then 


dx 


•  K 


dx 


at  P. 


(25b) 


In  this  casej  assume  that 


s(t) 


(t-t”r 

2t 


s(t")  + 


it-t 


5» 


im\5  _ 

8(t") 


+  ••• 


Q*(t)  - 


(t-t'”)  Q*(t'') 


+ 


2t 


(f»(t-) 


+ 


•  •  • 


t  a  t'"  (24  ) 


where  clearly 


i  (I*(t'") 
83(0  Xi(t’") 


(24  a) 


Differentiation  of  the  second  of  (I5)  along  s  gives 


6 


In  the  equations  which  follow,  all  quantities  must  be  evaluated  at  x  ■  s(t). 
The  similarity  of  this  portion  o£  the  proof  with  that  corresponding  one  of 
Theorem  IV  of  [1]  will  be  readily  noted. 
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whlle  the  first  of  (15)>  written  in  differential  form,  is 


a*  ®  dt 


dx 


s  + 


5i 

at 


0 


With  (7)  and  (ll),  Eqs.  (26a)  can  be  rewritten  as 


(26a) 


(26b) 


Writing  each  of  Eqs.  (25b)  for  solutions  I  and  2  and  subtracting  the  results 
in  each  case  one  obtains,  with  the  aid  of  (25b), 


a.2  ”  ax2 


0  at 


(26c) 


•J 

The  same  process,  applied  to  Eq.  (25)  now  gives'  at  Pg,  with  the  aid  of  Eqs. 

^  After  use  of  Eq,  (7)  to  give  ^  f  I  .  k  - ^  +  f(-^^.— ) 

axat  ■  ax^J  ^  3^3  ^  ^ 

and  similarly  for  Eq.  (ll)  and  _ L  . 

dxat 
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(23)  and  (26): 


liOll  .k  H  — _  «J*-d^)-0  tCd  -8  )  (27) 

^3  “Sn  *‘sm  ^3  “  ^'^2^1^  Psm*^®2  *1^ 


Id  view  of  Eq.  (23b)  and  (26c )>  the  argument  lnmedlately  following  Eq.  (22a) 
now  requires  that  the  left-hand  side  of  (2?)  be  non-positive;  this  again  leads 
to  a  contradiction  unless  both  equality  signs  hold  in  Eqs.  (24b).  The  case 

Q*  .  ;  5*  -  a*  I 

>  at  t  -  t'"  (28) 

Oj^  ■  Sg  ;  Sj^  a  ip  ;  •]_  “  *2  ) 


oust  therefore  still  be  investigated,  as  well  as  subsequent  special  cases  of 
this  type.  These  are  treated  by  means  of  further  differentiation  of  Eqs.  (23)« 
(26a),  (7)  and  (11),  along  s;  but  inspection  reveals  that,  in  general,  the 
result  will  be  (at  Pg): 


*'lai  2n+l  ^ 


dx' 


Sm 


) 


dx‘ 


2n+l 


dt 


P  i 


) 


dt 


IH-1 


(29a) 


ax“ 


&°*r  -T  ) 
S2  si^ 


a  0  ;  m  a  0,  1,  2, 


2n 


(29b) 


when 


d“(Q*-qp 

- a  0  ;  m  a  0,  1,  2,  •••,  n  -  1  (29c) 

dt 
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- «— ±- -  0  ;  m  -  0,  1,  2,  n  (29d) 

dt 

Choose  n  as  the  smallest  integer  for  which  at  least  one  term  of  the  right-hand 
side  of  (29a}  does  not  vanish^  and  then  note  that,  just  as  before,  this  leads 
to  a  contradiction  with  the  non-positive  character  of  the  left-hand  side  of 
(29a )•  The  proof  of  the  theorem  is  thus  complete. 

It  should  be  noted  that  the  converse  of  this  theorem  is  false,  that  is, 
the  validity  of  Eqs.  (1?)  does  not  necessarily  imply  the  validity  of  (16).  As 
a  corollary  of  this  theorem,  the  statement 


t 


12 


S  t 


11 


(50) 


follows  directly,  where  t^  is  (cf.  Eqs.  1^)  the  tisw  at  which  the  entire  slab 
has  melted. 

Another  type  of  bound  on  the  solution  will  now  be  established,  namely 
one  obtained  by  comparison  of  the  rate  of  melting  in  the  problem  defined  in 
Section  3  (in  which  the  sielted  portion  remains  stationary)  with  that  in  the 
companion  problem  in  which  the  melted  portion  is  instantaneously  removed.  This 
is  done  by  means  of  the  following: 

Theorem  VI  Consider  two  pairs  of  functions  Tj^(x,t),  Sj^(t)  and  Tg(x,t), 
such  that  the  pair  T^^,  s^  is  a  solution  of  Eqs.  (7)  to  (14)  with 
Q*(t)  a  0,  and  the  pair  Tg,  s^  satisfies  Eqs.  (7),  (8),  (9)«  (10)  and  (14)  in 
82(1)  <  X  <  1  as  well  as  the  following  two  equations: 

T2(s2(t),tl  -  T^ 

ar.  ds  (t) 

•‘‘Sm  dT  -  Psm^  dT’ 


(51) 
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where  the  same  function  Q(t}  Is  used  in  Eqs.  (lO)^  (12)  and  in  the  second  of 
(31);  then 


s^(t)  s  SgCt)  ,  ^  ^ 


(32a) 


and 


Tj^(x,t)  i  T2{x,t)  ,  t  i  t^  ,  82{t)  <  X  <  L  (52b) 


This  theorem  therefore  states  that,  under  the  same  heat  input  history,  a  more 
repld  advance  of  the  Interface  occurs  when  the  melted  material  is  instantane¬ 
ously  removed  (problem  2)  than  when  It  remains  stationary  (problem  l).  This 
conclusion  was  already  reached  in  [3I  the  special  case  of  Q  ■  constant;  it 
was  also  noted  there  that  for  very  short  times  after  the  onset  of  melting  the 
solutions  to  the  two  problems  are  identical. 

To  prwc  this  theorem,  we  start  with  the  heat-balance  equation  (16)  for 
problem  1  (with  Q*  ■  O)  and  the  corresponding  equation  for  problem  2,  namely 


t  L 

/  q(t)  dt  -  /  Hep  dx  +  (Pcm^  +  H  )  s  (t) 
0  s,(t)  Sm  m  2 


(53) 


Assume  now  that  s^  2  8^  (Fig.  ^);  then  subtraction  of  (33)  (16)  gives 


Sl(t)  L  ®l(t) 

/  IL  dx  +  /  (Hg.-H  )  dx  -  /  H  dx  +  (p  A+H  )[8j(t)-s  (t)]  .  0 
0  s^(t)  82(t)  Sm  m  1  2 


(33a) 


Now,  on  X  ■  Si(t),  T,  «  T„  and  T„  s  T  ,  so  that  T,-T„  2  0  (and  therefore 
'  l  im  2  m  12 

i  0)  for  the  entire  range  8j^(t)  s  x  ^  L,  by  Theorem  II,  Furthermore 
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fll(t) 

;  H  dx  s  H  [s,(t)-8  (t)]  (55b) 

B^(t)  m  1  2 

Therefore  Che  left-hand  side  of  (^5^)  is  composed  only  of  non-negative  terms 
and  as  a  consequence  this  equation  cannot  be  satisfied.  Hence  the  differ¬ 
ence  cannot  be  positive  (and  thtis  neither  can  theorem 

is  established. 

Problem  2  is  thus  seen  to  provide  an  upper  bound  to  the  solution  of 
problem  Ij  a  fortiori,  upper  bounds  to  the  solution  of  problem  2  are  also  upper 
bounds  to  the  solution  of  the  present  problem.  A  method  for  constructing  bounds 
for  problem  2  was  developed  in  [1]. 

Between  the  two  extremes  of  instantaneous  removal  and  of  stationary  melt 
there  may  be  defined  intermediate  problems  corresponding  to  finite  rates  of 
ablation.  It  may  be  conjectured,  as  an  extension  of  the  last  theorem,  that  a 
monotonic  relationship  exists,  in  the  problems,  between  the  rate  of  ablation 
and  the  rate  of  advance  of  the  solid-liquid  interface;  this  question  is  however 
not  examined  here. 

5.  Construction  of  Bounds;  Example 

The  procedure  for  the  use  of  the  bounds  previously  derived,  in  estimating 
the  solution  to  an  actual  problem,  consists  essentially  of  constructing  a  solu¬ 
tion  of  Eqs.  (7)  to  (14),  with  a  heat  input  ^J(t)  at  x  »  0  which  may  or  may  not 
equal  the  prescribed  input  Q(t),  disregarding  however  the  last  of  Eqs.  (15)« 

The  latter  equation  is  then  used  to  calculate  Q''^(t},  and  Theorem  V  Insures  that 
either  an  upper  or  a  lower  bound  has  been  found  in  a  range  0  <  t  <  t^^,  according 
to  whether  the  relations 
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Q*(t)  i  0  ;  Q(t)  +  Q*(t)  k  Q(t) 


(5*^) 


are  always  satisfied  or  always  violated  in  that  ranga«  It  will  be  noted  that, 
though  a  little  more  complicated  in  practice,  this  procedure  is  quite  analo¬ 
gous  to  that  described  in  [1]  for  the  problea  in  which  the  melt  is  instantane¬ 
ously  removed,  and  therefore  will  be  discussed  only  briefly  here. 

It  is  convenient  to  split  the  8olutl.on  in  two  parts,  the  first  pertain¬ 
ing  to  the  solid  portion  and  the  second  to  the  liquid  portion.  The  solution 
for  the  solid  is  easily  derived  by  consider'ing  s  fictitious  extension  to  a  slab 
of  the  original  thickness  L,  under  an  arbitrarily  chosen  fictitious  heat  input 
Q(t)  at  X  ■  0.  It  is  shown  in  [l,3l  that  the  relation  T„(s,t)  m  T  then  leads 
to  the  following  equation  for  8(t),  for  the  case  of  constant  properties  in  the 
solid: 


y  g(y-y^)  •  -t2nL'+5(y)]2/yj^  -[2(n+l)L’-5(y)]‘^/y 


y+1  f(y-y,  •  -[2nL'+?(y)]^/y,  -l2(n+l)L'-5(y)]P/y, 

-  2a  -  /  Lie  S  e  dy  (55) 

0  '^'1  n-0 


with  the  following  dimensionless  notation: 


(35») 


where  is  a  constant  reference  beat  input.  In  deriving  this  equation,  the 
convenient  choice  Q(t)  >  Q(t)  for  t  <  t^  «aas  made. 
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Tbe  solution  for  Che  liquid  portion  now  requires  the  determination  of 
Che  solution  of  Eq.  (ll)  In  Che  region  0  <  x  <  s{t)f  t  >  C^,  satisfying  the 
conditions 

Tl  \ 

^  .  Q(t) 

where  the  function  8(t}  Is  given  by  Eq.  (3^).  Since  s(t)  will  be  In  general 
known  only  numerically,  a  nunterlcal  solution  of  the  problem  for  the  liquid  Is 
probably  the  most  qpproprlate;  the  needed  solution  of  the  heat-condltlon  equa¬ 
tion  for  a  domain  bounded  by  boundaries  moving  In  a  prescribed  manner  can  In 
fact  be  conveniently  carried  out  by  such  methods.  For  the  purpose  of  Illus¬ 
trating  Che  bounds,  however,  analytical  expressions  valid  for  short  times  will 
be  derived  by  the  Integral-equation  method  of  [jlt  According  to  that  method, 
T^(x,t)  can  be  obtained  as  the  temperature  In  the  region  0  <  x  <  *,  under  the 
heat  input  Q(t)  at  x  -  0  and  Initially  (l.e.  at  t  ■  t^)  at  an  Initial  tempera¬ 
ture  distribution  T,  (x,t  )  -  X  6(x),  vThere  8(1)  ■  1  and  X  ■  x/^H.t.  The  first 
of  Eqs.  (^)  then  gives  the  following  Integral  equation  for  0(x),  for  the  case 
In  which  the  properties  of  the  liquid  axe  uniform  but  not  necessarily  the  same 
as  those  of  the  solid: 


t  s  t 


m 


(56) 


g 

°  This  Is  a  special  case  of  the  a ithod  of  [^];  sure  generally,  one  may  take 
the  region  0  <  x  <  (L^  2  L),  under  the  conditions  listed  above  and  In 

addition  under  a  heat  unpuc  at  x  >  to  be  suitably  specified. 
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with  the  following  additional  dlmenslonleaa  notation: 


Once  Eq.  (37)  ^  aolved,  the  last  of  Eqa.  (I3)  is  used  to  determine  Q*(y)> 
Consider  now^  for  simplicity,  the  problem  corresponding  to  Q(t)  ■  Q^, 
a  constant,  so  that  or  ■  1  and  so  that  the  exact  solution  of  [3],  for  short 


times,  is: 


5(y)  -  ^  y^^^  "l^y  *  0(y5'^) 

f(y)  -  -  f  ^  y  * 
e(x)  -  1  -  K  /«  X  +  o(x^) 


To  establish  the  bounds,  let  then 


ff  ^  2  1/2  m 

Hy)  “  ~  „  y  *  *y 


where  a  is  a  constant  to  be  discuss ed  presently.  Substitution  of  this  and  of 
f(y)  ■  1  into  Eq,  (35)  gives 


?(y)  •  ^  «y^^^  +  0(y) 
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so  that  clearly  for  short  times  this  is  an  upper  or  a  lower  bound  according  to 
whether  a  >  1  or  a  <  1,  the  exact  result  corresponding  to  a  ■  1.  The  short- 
time  solution  of  (37)  is,  with  T  b  1, 

e(x)  -  1  +  o(x)  (39b) 

and  finally  the  last  of  Eqs.  (I3)  gives 

J  (a-1)  +  0{y)  (39c) 

Clearly,  for  y  «  1,  Q*  b  0  If  a  >  1,  as  it  should  be  for  the  exact  solution; 
furthermore,  >  0  or  Q*  <0  according  to  whether  a>lora<l.  Since  here 
Q'(t)  a  Q(t),  the  former  of  these  possibilities  correspanda  to  an  upper  bound, 
and  the  latter  to  a  lower  bound,  according  to  relations  (3^);  this  is  plainly 
in  agreement  with  the  melting  rates  found  in  Eq.  (39^)* 

In  conclusion,  it  may  be  remarked  that  a  simple  bound  to  the  solution 
of  the  problem  defined  by  Eqs.  (7)  to  (14)  is  easily  found  in  the  special  case 
in  which  the  properties  of  the  solid  and  liquid  are  uniform  and  equal.  If 
again  Q(t)  b  and  for  a  semi-infinite  solid,  one  may  take  the  temperature 
both  before  and  after  melting  as 

T(x,t)  B  — -  ierfc  ,  0  <  x  <  •  (40) 

It  is  then  easy  to  show  that  the  relation  T  b  X  gives 

n 


1  -  /«(l+y)  ierfc  ( g//T+y) 


(40a) 


and  -chat 


where 


Qq  m  dy  Wl+y  erfc  (5//I+y) 


m  M 


ta 

2£ 


(40b) 


(41) 


This  can  be  readily  shown  to  correspond  to  an  upper  bound  (except  In  the  case^ 
here  trivial^  of  m  >  in  which  it  is  an  exact  volution  with  Q'*’  m  O):  for 
examplej  for  y  «  1,  Eqs.  (40a)  and  (40b)  give 


i 


1 

a/« 


2!kl 

% 


(42) 


in  agreement  with  the  conclusions  of  Theorem  V 
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